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1 . 0  INTRODUCTION 

The  towboat  manuevering  simulator  is  based  on  the  integration 
in  time  of  the  differential  equations  which  describe  the  motions 
of  the  towboat  and  barge  string  in  three  degrees  of  freedom,  i.e., 
yaw,  sway  and  surge.  The  theoretical  background  for  these  equations 
are  presented  in  Reference  1,  "The  Prediction  of  River  Tow  Maneu¬ 
vering  Performance,"  U.S.  Coast  Guard  Report  No.  CG-D-32-78.  This 
reference  presents  the  basic  equations  and  a  complete  set  of  hydro- 
dynamic  coefficients  for  a  representative  towboat  and  barge  train. 
These  coefficients  were  obtained  by  model  tests. 

This  section  of  the  simulator  documentation  provides  a  des¬ 
cription  of  the  basic  equations  of  motion  included  in  the  simulator 
and  the  relationships  used  to  determine  external  forces  and  moments 
due  to  a  bow  thruster  and  wind.  These  equations  are  completely 
general  in  nature  and  could  be  used,  with  the  proper  hydrodynamic 
coefficients,  to  describe  the  maneuvering  of  vessels  other  than 
towboats . 
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2.0  MATHEMATICAL  MODEL 

The  mathematical  model  for  the  maneuvering  of  a  river  tow 
consists  of  the  coupled  differential  equations  in  three  degrees 
of  freedom  (yaw,  sway  and  surge)  which  describe  the  motions  in 
the  X,  Y  plane  and  the  complete  set  of  hydrodynamic  coefficients 
and  external  forces  which  are  required  in  order  to  numerically 
integrate  these  equations .  There  are  also  auxiliary  equations 
which  describe  the  response  of  the  steering  and  propulsion  sys¬ 
tem  to  external  inputs. 

A  complete  set  of  three  coupled  differential  equations  with 
all  of  the  terms  necessary  to  simulate  normal  maneuvers  of  sur¬ 
face  ships  are  presented  in  Reference  2.  These  equations  have 
been  used  by  HYDRONAUTICS,  Incorporated  to  calculate  the  maneuver 
trajectories  for  a  wide  range  of  surface  ship  types  in  deep  and 
shallow  water.  These  equations  are  based  on  a  more  complete  set 
of  equations  developed  by  the  U.S.  Navy  for  the  simulation  of 
submarine  motions  in  six  degrees  of  freedom.  The  equations  in 
Reference  2  differ  from  other  sets  of  equations,  such  as  those 
of  Reference  3  used  to  describe  surface  ship  maneuvers  pri¬ 
marily  in  the  way  higher  order  terms  are  introduced.  The  equa¬ 
tions  of  Reference  3  use  a  Taylor  expansion  which  results  in 
odd  functions  being  represented  by  linear  and  cubic  terms.  The 
equations  of  Reference  2  are  a  square  absolute  representation  for 
higher  order  terms  so  that  odd  functions  are  represented  by 
linear  and  square  terms.  At  large  drift  angles,  which  is  a  likely 
operating  condition  for  a  river  tow,  the  forces  and  moments  are 
dominated  by  cross  flow  drag  which  is  proportional  to  velocity 
squared.  Thus,  it  is  better  to  use  equations  in  which  forces  and 
moments  are  proportional  to  velocity  squared  rather  than  cubed. 

As  a  result  ,  the  equations  presented  in  Reference  2  were  selected 
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to  form  the  basis  for  the  equations  which  describe  the  maneuvering 
of  a  river  tow.  There  are  a  number  of  modifications  which  must 
be  made  to  these  equations.  These  modifications  and  the  resulting 
set  of  equations  are  described  in  the  following  paragraphs. 

The  following  equations  in  three  degrees  of  freedom  are  re¬ 
ferred  to  a  right-hand  prthogonal  system  of  moving  axes,  fixed  in 
the  body,  with  its  origin  normally  located  at  the  center  of  mass 
of  the  body.  The  positive  direction  of  the  axes,  angles,  linear 
velocity  components,  angular  velocity  components,  forces  and  moments 
are  given  in  Figure  1.  Unless  otherwise  indicated  in  the  Notation, 
the  numerical  values  for  the  hydrodynamic  coefficients  used  with 
the  equations  are  for  the  ship  propulsion  point  (n  *  1.0).  The 
equations  are  written  in  terms  of  the  complete  barge  flotillia 
towboat  configuration.  Thus  the  values  of  the  coefficients  embrace 
the  interaction  effects  between  rudder  and  hull,  propeller  and  hull, 
and  propeller  and  rudder  as.  determined  from  the  model  tests  of  the 
complete  configuration. 

An  important  consideration  in  the  maneuvering  of  a  river  tow 
is  the  effect  of  current  which  can  vary  significantly  along  the 
length  of  the  tow.  As  a  result,  it  is  necessary  to  introduce  the 
effect  of  the  current  velocity  into  the  mathematical  model.  The 
approach  adopted  was  to  define  the  hydrodynamic  terms  in  the  equa¬ 
tions  based  on  the  relative  velocities  and  yaw  rate  between  the 
hull  and  the  fluid  rather  than  the  inertial  velocities  and  yaw 
rate.  The  relative  velocities  and  relative  yaw  rate  can  be  cal¬ 
culated  by  the  vector  addition  of  the  Inertial  velocity  and  iner¬ 
tial  yaw  rate  and  the  current  velocities  and  current  yaw  rate.  In 
the  numerical  integration  the  procedure  is  to  define  a  matrix 
of  current  speeds  (Ucij)  and  directions  (?cl ^ )  at  points  on  the 
X,  Y  plane.  Based  on  the  location  of  the  bow  (XB,  Y^)  midships 
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(X  .  Y)  and  stem  (X  .  Y„)  of  the  tow,  an  interpolation  in  the 
current  speed  and  direction  matrix  is  carried  out  to  obtain  the 
current  speed  and  direction  at  the  bow  (UCB ,  i|>CB) ,  midships 
(UCM,  and  stem  (Ucg,  4ics> .  The  following  relationships  then 

apply: 


UCB 

■ 

UCB 

cos  (i|>CB-t|>) 

VCB 

-°CB 

sin 

UCM 

- 

UCM 

COS  (fgg-f) 

VCM 

*  UCM 

sin 

ucs 

- 

ucs 

COS 

• 

vcs 

-  BCS 

sin 

^cs*^ 

if 

UCB 

+  “cm-*-  ucs 

+  v_,  +  v„„ 
CM  CS 

uc 

3 

V 

c 

3 

r 


c 


i 


(l) 


“r 


u-u. 


'R 


v-v„ 


R 


r-r. 


In  this  procedure  the  mean  longitudinal  and  lateral  current  velo¬ 
city  in  the  body  axis  system  is  obtained  from  the  average  of  the 
values  at  the  bow,  midships  and  stem.  The  variation  of  the 
lateral  velocity  along  the  tow  is  accounted  for  by  the  apparent 
current  yaw  rate  defined  by  the  difference  in  the  lateral  velo¬ 
cities  at  the  bow  and  stem  divided  by  length.  This  assumes  the 
lateral  velocity  varies  linearly  from  bow  to  stem.  If  this  is 
not  the  case  a  more  complex  relationship  would  have  to  be  intro¬ 
duced  . 

The  equations  of  motion  formulated  for  a  river  tow  are  as 
follows: 


n 
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■  \  cos 


0 


YAWING  MOMENT 


Propulsion  Ratio  -  (n/u_)/(n  /u 
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The  foregoing  equations  of  motion,  as  noted  previously,  are  pat¬ 
terned  after  the  quasi-steady  state  equations  of  Reference  2.  The 
non-dimensional  hydrodynamic  coefficients  which  comprise  the  basic 
equations  are  considered  to  be  independent  of  speed  (Froude  number) . 
This  assumption  is  valid  since  river  tows  always  operate  at  low 
Froude  number. 

These  equations  differ  from  the  equations  of  Reference  2  in 
the  following  details: 

a.  They  are  written  in  terms  of  the  relative  velocities 
and  yaw  rate  to  allow  the  introduction  of  varying  current  as 
discussed  above. 

b.  Terms  are  included  for  steering  rudders  (steering  rudder 
angle  »  5g)  and  flanking  rudder  (flanking  rudder  angle  -  6p) . 

c.  Terms  are  included  for  twin  propellers  which  nay  operate 
at  different  RPM's  and  different  directions  of  rotation.  The 
turning  moment  due  to  differential  thrust  is  included  in  the  yawing 
moment.  Equation  (4) . 

d.  The  forces  and  moments  generated  by  the  rudders  are  based 
on  a  velocity  defined  by 

Up*  *  <dun*  +  eDnuR  +  fD2n2) 

which  is  a  function  of  the  relative  axial  velocity,  uR,  and  pro¬ 
peller  RFM,  n.  The  constants  d,  e  and  f  depend  on  the  sign  of 
tig  and  n.  This  allows  a  proper  representation  of  the  rudder 
forces  and  moments  at  zero  speed  and  finite  propeller  RFM. 

In  realistic  maneuvers,  river  tows  operate  both  ahead  and 
astern  and  in  some  cases  at  large  drift  angles.  In  order  to  pro¬ 
perly  represent  the  hydrodynamic  forces  and  moments  which  act  in 
such  conditions,  different  sets  of  hydrodynamic  coefficients  are 
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used  depending  on  the  relative  drift  angle  0R.  As  will  be 
noted  in  Reference  1  in  which  the  hydrodynamic  coefficients 
from  the  model  tests  are  presented,  most  coefficients  depend  on 
the  direction  of  motion,  i.e.,  ahead  270°  *  c  90°  or  astern 
90°  *  8r  *  270° .  In  order  to  obtain  a  better  representation  of 
the  steady  sideforce  and  yaw  moment  at  drift  angles  near  90  and 
270  degrees  certain  coefficients  have  an  additional  value  when 
30°  -  -  150°  or  210°  -  0R  5  330°. 

As  noted  in  the  introduction  to  this  section,  a  complete  set 
of  hydrodynamic  coefficients  for  Equations  2,  3  and  4  for  a  tow- 
boat  and  barge  string  are  presented  in  Reference  1.  The  simulator 
as  presently  configured  has  the  hydrodynamic  coefficients  as  con¬ 
stants  independent  of  water  depth.  Thus,  if  operations  in  shallow 
water  are  to  be  simulated,  hydrodynamic  coefficients  applicable  to 
the  appropriate  depth  should  be  used.  As  a  general  rule,  the  water 
can  be  considered  as  deep  if  it  exceeds  the  draft  of  the  tow  by 
a  factor  of  2.5  or  3. 

The  equations  of  motion  (i.e. ,  equations  2,  3  and  4)  are 
solved  stepwise  in  time  in  the  computer  program.  In  the  program 
the  time  step,  DT  is  set  equal  to  1.0  seconds.  For  a  long,  re¬ 
latively  slow-moving  river  tow,  a  longer  time  step  (2  to  4  seconds) 
could  be  used  without  significantly  affecting  the  results  of  most 
maneuvers.  The  1.0  second  time  step  was  chosen  more  on  the  basis 
of  the  updating  rate  for  the  control  display.  It  was  considered 
desirable  to  minimize  the  time  delay  between  a  control  input  and 
the  display  response.  A  1.0  second  time  step  was  considered  rea¬ 
sonable  on  this  basis. 

At  each  time  step  in  the  solution  of  the  equations  of  motion, 
Che  current  velocity  at  the  low,  midship  and  stem  of  the  tow  is 
determined.  This  calculation  is  carried  out  by  subroutine  CURT. 
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In  the  Input  to  the  program,  the  current  speed  and  direction  is 
specified  at  up  to  30  lines  or  stations  across  the  river.  These 
may  be  arbitrarily  spaced  along  the  river.  The  river  bank  is 
defined  in  the  visual  display  by  lines  joining  the  end  points  of 
consecutive  stations.  As  a  result,  stations  should  be  concentrated 
in  the  area  of  bends  or  rapid  changes  in  river  width.  The  current 
speed  and  direction  are  also  likely  to  change  rapidly  in  such 
areas. 

At  each  station,  the  current  speed  and  direction  is  specified 
at  8  evenly  spaced  points . (see  Figure  13  of  Volume  I) .  Subroutine 
CURT  determines  the  current  speed  and  direction  at  bow,  midship 
and  stern  by  searching  the  stations  in  sequence.  At  each,  station, 
the  smallest  distance  between  the  vessel  point  and  one  of  the  8 
station  points  is  determined.  Initially  this  process  is  repeated 
until  the  minimum  distance  point  is  found.  The  current  defined  at 
the  closest  point  is  then  used  in  calculations.  After  the  initial 
time  step  the  search  is  carried  out  locally  around  the  station 
point  found  to  have  the  minimum  distance  during  the  previous  time 
step.  This  local  search  technique  saves  significant  computer  time.  . 


t 

i 


Jt 
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3.0  EXTERNAL  FORCES 


3.1  Bow  Thruster 

The  towboat  maneuvering  simulator  contains  the  provision 
for  a  bow  thruster.  The  forces  and  moments  generated  by  the  bow 
thruster  are  represented  by  the  following  relationships : 

Thruster  Lateral  Force  - 

FYo  *  *  (Percent  Output)  *  THMAX  *  TRPM 

FYq  •  Y  force  at  zero  forward  speed 

Percent  output  ranges  from  -1.0  to  1.0 

The  lateral  forces  of  a  tunnel  type  bow  thruster  decreases  with 
forward  speed.  Thus,  the  lateral  force  introduced  into  the  lateral 
force  Equation  (3)  is  given  by: 


where 


FXP(3)  -  Fy  -  Fy  (TRC(l)  +  TRC(2)  *  Cos  (2r  u/Vj))  (5) 

for  u/Vj  *  0.5 

-  Fy  <(TRC(1)  -  TRC(2))  +  TRC(3)  *  (U/Vj-0.5)) 

for  u/Vj  >  0.5 


axial  speed  of  vessel 

thruster  jet  or  exit  velocity 

<Ft  /PA)* 
o 

mean  density  of  water 


A  ■  thruster  Disc  Area 
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Typical  constants  for  a  tunnel  bow  thruster  based  on  Reference  S 
would  be: 

THMAX  -  10,000  lbs 
TRC(l)  -  0.65 
TRC(3)  -  0.18 

CTHV  -  1/pA  -  0.025 
Thruster  Yawing  Moment 

Because  of  interaction  forces  with  the  hull  the  yawing  moment 
due  to  the  thruster  does  not  change  with  forward  speed  in  the  same 
way  as  the  lateral  force.  The  yawing  moment  due  to  the  thruster 
which  is  introduced  into  the  yawing  moment  equation  is  given  by: 

N  -  Fy  Xt  (1  -  (1  -  0.67  u/VJ) ^1  -  yl-j)  (6) 

N  -  Yawing  Moment  due  to  thruster 
X£  -  Longitudinal  location  of  thruster  -  YCSP(3) 


3.2  Wind  Forces  and  Moments 


The  forces  and  moments  acting  on  the  towboat  and  barges  due 
to  wind  were  calculated  using  methods  and  data  given  in  Reference 
4.  The  forces  and  moment  at  a  given  time  are  calculated  using  the 
instantaneous  relative  wind  velocity. 


The  axial  and  lateral  forces  and  yawing  moment  due  to  wind 
are  given  by: 


X  -  -ip  C  V  2A„  cos2  0  -  FX(2) 

w  xw  we  px  we  v 

Y„  •  ipC  V  2A  sin2  0  -  FY(2) 

w  iKa  yv  we  py  we 


(7) 

(3) 


N 


“  Vic, 


FN(2) 


(9) 
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where 


C  and  C 
xv  yw 

A  and  A 


X  and  Y  are  longitudinal  and  lateral  forces  (in  ship 

w  w  coordinates)  due  to  wind 

Nw  is  yawing  moment  about  the  c.g.  due  to  wind 

Pa  is  air  mass  density 

C  and  C  are  longitudinal  and  lateral  drag  coefficients 
xv  yv 

A  and  A  are  longitudinal  and  lateral  projections  of 
px  py  above  water  hull  and  deckhouses 

Vwe  is  effective  wind  velocity  *  VwR 

B  is  angle  measured  from  the  ship  x-axis  to  the 

we  effective  wind  velocity  vector  *  BW 

x-  is  longitudinal  position  of  the  center  of 

cp  pressure  of  the  wind  force 

The  effective  wind  velocity,  including  the  effect  of  the  boundary 
layer  at  the  water  surface,  deduced  empirically  in  Reference  4. 
is: 


(xrfr)  <V„I+V> 

/  z  \0'15  r 

Writ]  [<«  +  v 


cos  8  )2  +  (v  +  V 
w  xr  N 


sin  8  )2 
w  w 


where 


z  is  vertical  center  of  area  of  the  ship 

ca 

V  and  V  are  longitudinal  and  lateral  components  of  the 
7  effective  wind  velocity 


u  and  v 


are  components  of  ship  velocity 
is  actual  wind  velocity 

is  wind  heading  angle  measured  from  ship  x-axis 
to  actual  wind  velocity  vector 


1 
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The  angle  8W#  Is  given  by 


The  center  of  pressure  is  determined  from  the  longitudinal  center 
of  area,  using  empirical  methods  and  data  from  Reference  4  and 
is: 


Xlcp  "  LBP 


-  (0.001  +  0.01  AR)  (Is  ’I  -  90) 
UBP  w 


where 


JBP 


ca 


is  ship  LBP 

is  longitudinal  position  of  the  center  of  lateral 
area  *  XLCA 


S  ' 
w 

AR 


is  the  value  of  8wa  in  degrees  -  BW 
is  hull  aspect  ratio,  AR  -2Apy/LBp2 


AS 


Based  on  the  data  given  in  Reference  4  and  estimated  above 
water  hull  and  deckhouse  shapes,  the  following  wind  force  co¬ 
efficients  would  apply  to  the  towboat  and  barges  configuration 
for  which  hydrodynamic  data  are  reported  in  Reference  1. 


lbp 

fli 

745 

ft 

Xc. 

AR 

m 

m 

58.8 

0.0167 

ft 

C 

XV 

m 

1.0 

G 

yv 

m 

1.0 

p* 

m 

1520 

ft* 

V 

m 

4635 

ft* 

2ca 

m 

6.2 

ft 

. . M 
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3.3  Other  External  Forces 

At  this  time  no  other  external  forces  are  included  in  the 
simulator.  In  the  future,  othev  external  forces  such  as  bank 
suction,  interaction  with  passing  vessels  and  the  effects  of  moor~ 
ing  lines  could  be  included. 
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4.0  CONTROL  SYSTEMS 

The  control  systems  included  in  the  towboat  maneuvering 
simulator  include  steering  control  for  the  steering  and  flanking 
rudders,  RPM  control  for  the  port  and  starboard  propellers  and  the 
bow  thruster  output  control.  The  command  inputs  for  all  of  these 
controls  are  provided  thru  the  computer  A  to  D  input  from  the  con¬ 
trol  station.  In  the  simulator,  the  response  to  the  control  inputs 
is  at  a  constant  rate.  The  rates  used  based  on  Reference  6  are 
as  follows: 

Steering  and  Flanking  Rudders  5  deg/ sec 
Propeller  RPM  20  RPM/ sec 

Thruster  Output  10  percent  thrust  change/ sec 

The  response  to  the  control  system  inputs  is  considered  re¬ 
presentative  of  a  typical  towboat.  If  required  in  the  future,  a 
more  complex  representation  of  the  response  to  control  Inputs  could 
be  modeled  in  the  simulator. 
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